Time drift of cosmological redshifts and its variance 
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The contribution of cosmological perturbations to the time drift of the cosmological redshift is 
derived. It is shown that the dominant correction arises from the local acceleration of both the 
emitter and the observer. The amplitude of this effect is estimated to be of the order of 1% of the 
drift signal at z — 2 — 4, but can easily be lowered down to 0.1% by using many absorption lines 
and quasars. 
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I. INTRODUCTION 

The increasing observational evidences that the expan- 
sion of the universe is accelerating (see Ref. 0, Q and 
reference therein) has stimulated a rising interest to the 
reconstruction of its expansion history. An important 
outcome of these theoretical studies is to clarify the sen- 
sitivity of observational tests to dark energy properties 
and to assess how each could be corrupted by extra-noise 
from other cosmological effects. Comprehensive investi- 
gations of these nuisances have been carried out on "stan- 
dard" tests, like SNIa, weak lensing, BAO, CMB, ISW 
or clusters of galaxies. 

In contrast, very few focussed on the time drift ef- 
fect that changes the observed redshift of an object as 
function of time. The recent claim that it may drive the 
conceptual design of instrument for next generation giant 
telescopes raised the need that similar attention should 
be paid to the theoretical ground of this novel technique. 
Interestingly, such an observation may lead to a better 
understanding of the physical origin of the recent accel- 
eration 0, Q| and to a determination of the dark energy 
equation of states [H as well as constraints on dark en- 
ergy models M or tests of the variation of fundamental 
constants 0,Q. 

As first pointed out by Sandage j|, in a homogeneous 
and isotropic spacetime, the time drift of the observed 



redshift is directly related to the Hubble function by 

z = (1 + z)Ho-H(z)=z(t) ,z) . (1) 

Given the most likely ranges of cosmological parameter 
values derived from observations, in a ACDM model the 
typical amplitude of the redshift drift is of order 5z ~ 
—4 x 1CP 10 on a time scale of St = 10 yr, for a source 
at redshift z = 4. This corresponds to a spectroscopic 
tiny velocity shift, Sv = cSz/(l + z), of Sv ~ 2.5 cm/s. 
Fig. [1] (left panel) shows the time drift as function of 
redshift for the standard ACDM model and a dark energy 
models with an equation of state changed by only 10% 
(w = —0.9), all other parameters being kept constant. 
Both curves have similar shape but the difference the 
drifts between a standard ACDM model and cosmological 
models tends to zero at hight redshift. Fig. [1] (right panel) 
depicts this difference for two models with either w = 
-0.95 or w = -0.98. 

The feasibility of this measurement is most challeng- 
ing and impossible with present-day astronomical facili- 
ties. However, it was recently revisited [lfj in the con- 
text of the new generation of Extremely Large Tele- 
scopes 1 (ELT), arguing that with such outstanding col- 
lecting areas one could measure velocity shifts of order 
Sv ~ 1 — 10 cm/s over a 10 year period from the ob- 
servation of the Lyman-alpha forest on QSO absorption 
spectra. In particular, it is one of the main science driver 
to design the COsmic Dynamics Experiment (CODEX) 
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spectrograph [Tl|, [l2| for the future European ELT (E- 
ELT). 

The performances of CODEX and its capability to 
measure a time drift of very distant objects were esti- 
mated using Monte-Carlo simulations of quasar absorp- 
tion spectra. The expected velocity accuracy of this ex- 
periment can be written as follows (see Ref. fill ]) 
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provided the absorption lines are resolved. S/N denotes 
the signal-to-noise ratio, for a pixel scale of 0.0125 A and 
A^Qgo is the number of quasars. Thus, spectroscopic mea- 
surements of about 40 quasars with S/N ~ 2000 ten 
years apart can reach a 1.5 cm/s accuracy. This is within 
the reach of a CODEX instrument mounted on a 60 — 80 
meter ELT by observing a 16.5 th magnitude QSO during 
2000 hrs 03. 

Many systematic effects that may spoil the time drift 
signal, such as Earth rotation, proper motion of the 
source, relativistic corrections etc., are discussed in 
Ref. [ll|. The acceleration of the Sun in the Galaxy 
seems more a serious problem because its amplitude may 
be of the same order than the cosmic signal. However, 
it has not been measured yet, so its nuisance is still un- 
known. On the other hand, subtle contaminations like ac- 
celerations produced by large scale structures have never 
been estimated in the error budget. The purpose of this 
work is to address this issue and to estimate whether it 
may hamper the cosmological interpretation of the time 
drift. 



direction of observation. O and E refer to the observer 
and emitter respectively (see the Appendix more precise 
definitions of all the variables involved in this equation). 

The first term at the right hand side of Eq. © clearly 
arises from the local position of the observer. The second 
term of Eq. ^ encodes Doppler effect due to the rela- 
tive motion of the observer and the source as well as the 
equivalent of the integrated Sachs- Wolfe term [3, EH - 
Eq. © is the analog of the (direction dependent) tem- 
perature anisotropy of the cosmic microwave background 
(CMB) compared to the mean CMB temperature. 



III. ESTIMATE OF THE VARIANCE 

The variance of £(e) can be split into contributions 
coming from the time dependence of the gravitational 



potentials, = (1 
acceleration, Q, = (1 



and from the peculiar 

z)[e.*]l. 

The estimation of C$ demands a full description of the 
time evolution of the potential, both at emission and ob- 
serving times. In the following we derive it and discuss 
its properties using the linear cosmological perturbation 
theory. The validity of this approach will be more thor- 
oughly addressed in the next section. 

Using the linear theory of structure growth, the density 
contrast can be split as S = D(t)e(x), where e(x) com- 
prises all details on the initial conditions. The growth 
rate D + is the growing solution of the equation 



D(t) + 2HD(t) = -H 2 n m (t)D(t), 



(4) 



II. COSMOLOGICAL PERTURBATIONS 

Eq. ([1]) relates the time drift of the observed redshift 
to the Hubble function, assuming a perfectly homoge- 
neous and isotropic Friedmann-Lemaitre spacetime. In 
the real universe, however, velocity terms arising from 
cosmological perturbations add up as noise contributions 
and increase the scatter of the redshift drift around its 
mean value. 

The distribution of the redshift drift can be predicted 
using the expression of z to first order in the cosmological 
perturbations. It is derived in the Appendix A of this 
work. At first order in the metric perturbations and in 
v/c it writes 



z = z(rj Q ,z) + ((x o ,r) ,e;z) 



with 

((x Q ,T] Q ,e;z) 



-$ z{r)o,z) + (l + z) 



e.v - * 



(2) 



.(3) 



This formula involves both Bardeen potentials, $ and VE", 
and the peculiar acceleration, v. A dot denotes a deriva- 
tive with respect to observer proper time and e is the 



where ^ m (< : ) is the time dependent reduced density pa- 
rameter for the gravitating matter (see ref. [l6| for de- 
tails). On sub-Hubble scales, Einstein equations imply 
that f = $ and A$ = %H 2 Vt m a 2 5. 

As the redshift increases, the dynamics of the universe 
is closer and closer to the one an Einstein-de Sitter Uni- 
verse. $ is therefore almost constant and $ is expected 
to vanish. This is no longer true at low redshift, when the 
cosmological constant (or the spatial curvature) starts to 
dominate. Instead, the time evolution of the potential 
writes 



* = H* {f(t) - 1] 



(5) 



where fit) = d In D + /d In a comprises the intrinsic evo- 
lution of the potential produced by the growing pertur- 
bations. In a flat ACDM, / is explicitly given by 



/(*)=!- 
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(6) 



where a = -f/oV^Ao and where 2-F1 is a hypergeometric 
function. 

Using Eq. §5§ it is then easy to express the r.m.s. of $ 
from the r.m.s. of the mass density fluctuations, as, as 
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derived from the Poisson equation. More precisely 



crs = 



d 3 k 

(2^)3 



Ps(k) 



1/2 



(7) 



where Ps is the power spectrum of the density contrast, 
(SkSk') = Ps(k)Sri{k + k'), and 6k are the Fourier modes 
of S. To estimate Ps, we adopt the prescription by Bond 
et al. [13 for the transfer function and the normalization 
<Tg = 1. The redshift dependence of the power spectrum 
is then the one of the growing mode, D + (z), normalized 
to unity at z = 0. 

Turning to the gravitational potential, it appears that, 
in the standard model of cosmology with a primordial 
spectrum of index n s ~ 0.95, the amplitude of the po- 



tential fluctuations is IR divergent. However, since the 
previous calculation is only valid for sub-horizon modes, 
it is necessary to introduce a cut-off for modes typi- 
cally beyond the Hubble scale. The expected potential 
fluctuations then drop to more realistic amplitudes of 
ct$ ~ 5 x 10~ 5 . It follows that, for a source at redshift 
z, the r.m.s. of i induced by the time variation of the 
gravitational potential is 



(ci 



[/(0) - 1] a $ 



(8) 



which is of order Ci ~ (1 ' 
indeed. 



x 10 5 i/n, a small number 




FIG. 1: (left) The time drift of the redshift as a function of the redshift of the source obtained from Eq. ([T]) for a ACDM model 
(solid line) and a model with a constant equation of state w = —0.9 for the dark energy (dashed line), (right) Amplitude of 
the r.m.s. of the systematic errors due to cosmic acceleration effects. The contribution of £ D (dashed line) is subdominant 
compared to the one of £ E (dotted line). The solid lines represents the difference between a standard ACDM model and 
cosmological models with either w = —0.95 (upper solid line) or w = —0.98 (lower solid line). 



The contribution of the peculiar acceleration, is less 
obvious to derive because we do not have a complete the- 
ory that describes the expected distribution of the local 
line-of-sight acceleration. However, in the cosmological 
linear theory not only are the metric components sup- 
posed to be small (as explicitly used above) , but also the 
density contrast and the velocity gradients (compared re- 
spectively to unity and H), see Ref. 16]. 

The Lyman-alpha forest is believed to be dominated 
by low density clouds of intergalactic medium, with in- 
dividual accelerations primarily triggered by large-scale 
structures. Assuming then linear theory holds in our con- 
text, the local acceleration writes i>i = —Hvi — di^/a, so 
that 



Cv(e,z) = {l + z) e l 



a 



(9) 



In terms of the dimensionless divergence 0(x) = diVi/aH, 
the linear continuity equation reduces to 9(t, x) = 



—f(t) 5(x) at linear order. This implies that the Fourier 
components of the velocity, density contrast and po- 
tential are related by k 2 Hvi(k) = f(t)aH 2 kiSk and 
fc 2 $ ! i(fc)/a = 3f2 m a H 2 ki5k/2. Using our previous es- 
timate of Ps, one easily derives the r.m.s. of the two 
contributions to Cv, 



2\l/2 _ 



= (1 



2^„o - /(0) 



Ho & (10) 



that depends on the emission time only through the fac- 
tor (1 + z), and 



where a 



2\l/2 



2 

d 3 fc 1 



H 2 (t)D+(t) a (11) 



;P(k,Z = 0). 



(2tt) 3 3fc 2 

These two terms are independent and should be 
summed quadratically. The resulting r.m.s. of de- 
picted on Fig. Q] (right panel) shows C E is the dominant 
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contribution at all redshifts. It rises to a percent level 
from z = to z = 4. At z ~ 4, Q ~ 0.5%, while &(e, z) 
is ten times smaller. Both terms have similar behaviour 
and are basically unchanged for any realistic flat cosmol- 
ogy having and effective w close to — 1 , but note that this 
is a priori not the case for any model. 



IV. DISCUSSION 

Assuming the cosmological time drift derived from 
QSO absorption lines by the Lyman-alpha forest may be 
contaminated by extra-acceleration of clouds by massive 
structures, it is legitimate to question the validity of the 
linear regime approximation used throughout this work. 

Let us first consider the acceleration of an absorbing 
Lyman-alpha cloud. On large scales, clouds are located 
inside filaments infalling towards massive clusters or 
super-clusters of galaxies. Assume, then, the acceleration 
is due to the gravitational attraction of a super-cluster 
with typical mass of order 10 15 M Q , localized at 10 Mpc 
from the cloud. The Newtonian acceleration is about 
<2N ~ 1-45 x 10 _15 km/s 2 . In comparison, the Hubble 
acceleration cHq is an ~ 6.8 x 10 _13 km/s 2 so that 
on/oh ~ 2 x 10~ 3 . This ratio may change by one 
order of magnitude, depending on the mass and length 
scales one may consider for clusters, super-clusters or 
filaments, but is always sufficiently small to keep the 
linear approximation valid. It is also worth noticing its 
amplitude is close to theoretical expectations derived 
in the previous Section. We therefore speculate the 
simple interpretation of our theoretical estimate as 
being primarily due to the acceleration of the nearest 
rich cluster is pertinent 2 . To confirm this and get more 
sophisticated description of accelerations the use of 
numerical simulations is indeed necessary. 

In practice, a time drift is not measured from a sin- 
gle absorption line but by averaging several lines spread 
over a spectral range AA defined by the spectrograph. 
If the acceleration of Lyman-alpha clouds is primarily 
driven by clusters of galaxies located around their neigh- 
borhood, then clouds are not dynamically independent 
and accelerations of closeby clouds are correlated. We 
are thus interested in the variance of i, averaged over a 
bound comoving distance A\ along the line of sight, 

z = / )°Y ■ 



2 Liske et al. (in preparation) also estimated the contamination 
of the drift signal produced by peculiar motions. In contrast 
to our analysis done in a full General Relativity context, they 
simply used Special Relativity formalism. Note that FlU derived 
the peculiar acceleration of strong gravitational potentials like 
clusters of galaxies but to predict the peculiar velocity drift over 
several decades produced by nearby systems on a test particle. 
Both results agree with our predictions. 



It is related to the variance from correlations obtained 
on a single line by 

(z 2 )=a 2 (z,Az)(( 2 (z)) , 

where the coefficient a{z, Az) depends on the physical 
size over which the average is performed. Az is the red- 
shift range explored by the spectrograph at the mean 
redshift z: Az = (1 + z)AX/X. For a ACDM uni- 
verse, it corresponds to a comoving distance of Ax = 
D Ho [fi m0 (l + z) 3 + fiAo] _1/2 , with D Ho = 3000/1" 1 Mpc. 
a can be computed from the correlation of the accelera- 
tion field, 



(a(xi)a(x2)) = 



d 3 k 



3 iMxi-X2) 



P(k) 



as 



a 2 = 



(2tt) 3 3fc 2 
1 f dk z d 2 k± sinfc z Ax 



3fc 2 



k z A X 



P(k) , (12) 



where A% is the size of the comoving radial distance over 
which the average is performed. 




AA [A] 

FIG. 2: The coefficient a that enters in Eq. (I13|l as a function 
of the width of wavelengths, AA, on which the observations 
are average for several source redshifts, z = 4 (solid line), z — 
3 (long dashed line), z — 2 (dashed line) and z = 1 (dotted 
line) for ACDM with h = 0.7, f2 m0 = 0.3 and fi A o = 0.7. 

If we could naively split the Lyman-alpha forest along 
a line if sight into radial bunches of physically decoupled 
cloud systems, without correlated accelerations, 1/a 2 
would provide an estimate of the number of bunches. 
From an observational point of view, 1 /a 2 expresses the 
effective number of absorption line systems without cor- 
related acceleration probed by a spectrograph covering a 
wavelength range AA around the mean redshift z. It in- 
creases when the spectral coverage of the spectrograph 
increases (see Fig. [5]). For example, if z = 4, and 
AA = 100 A, 200 A, 500 A, that is Az = 0.1,0.2,0.5 
(see Fig. [2]) respectively, then a is 0.69, 0.55, and 0.38, 
and 1/a 2 = 2.1, 3.3 and 6.9. However, 1/a 2 only takes 
into account one light of sight. If we average over A^qso 
randomly selected lines of sight, then we expect that 



a(AX,z)N-^ 2 C(z) 



(13) 
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Hence, the spectral range of the spectrograph together 
with the number of lines of sight can easily drop the 
contribution of cosmological perturbation to the variance 
budget below a 0.1% level. 

It is interesting to notice that the theoretical values of 
1 /a 2 derived in the previous paragraph can easily be in- 
terpreted and predicted from simple physical arguments. 
The spectral range of a CODEX-like spectrograph as de- 
scribed in Ref. [u|. is AA ~ 500A. At a redshift of z — 4, 
it corresponds to Az ~ (1 + z)AA/5000A ~ 0.5 and to 
a comoving length of Ax ~ 300 Mpc. If we assume that 
the coherence scale of the velocity field is the typical size 
of the super-cluster (~30 Mpc), then CODEX can probe 
about 10 independent systems per line of sight. This is of 
the same order as l/a 2 = 6.9 for AA = 500 A discussed 
above, which confirms its interpretation as an effective 
number of uncorrelated cloud systems. It also simply 
explains why az cx a, which is nothing but the inverse 
square root of this number. 

V. CONCLUSION 

In order to measure the cosmological time drift of the 
redshift, many systematic effects will have to be under- 
stood. Besides the systematic errors of astrophysical ori- 
gin that may affect the observation of the Lyman-alpha 
forest, large scale structures will induce a dispersion of 
z. This work addresses this issue. 

First, we have derived the expression of the time drift 
of cosmological redshift at first order in the perturbation. 
This was then used to estimate its variance and then to 
demonstrate it depends on two main effects, the acceler- 
ations and the local gravitational potential at both the 
source and observer positions. 

The contributions at the observer position have not 



been discussed further. High precision astrometric obser- 
vations with GAIA will soon provide exquisite knowledge 
of the motion of the Earth in the Milky Way. It will pin- 
point all local acceleration terms with enough accuracy 
to remove this contribution easily fl9| . 

In contrast, the contributions at the source position 
are much more difficult to subtract. In the linear regime, 
we have shown that the gravitational potential contribu- 
tion is negligible while the acceleration of the source is 
typically of the order of 1% at z = 2 — 4. We argue a 
dominant contribution to this term is the acceleration of 
galaxy clusters near the source. 

In order to understand whether the amplitude of this 
variance can be reduced, we have estimated the effect of 
averaging the signal over several absorption lines. One 
can either profit from the total spectral range covered 
by the spectrograph to measure the drift from all lines 
detected along a line of sight, provided correlated accel- 
eration contributions are taken into account, or use the 
mean drift over many randomly selected lines of sight. 
The first option reduces the variance by the square root 
of the number of uncorrelated clouds systems along a line 
of sight, the second by the square root of the number of 
independent lines of sight. For an instrument having the 
current specifications of the CODEX spectrograph, it is 
then easy to drop the contribution of large-scale struc- 
tures to the total variance of the time drift down to a 
0.1% level. 
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APPENDIX A: TIME DRIFT AT FIRST ORDER 
IN THE PERTURBATION 

The redshift is defined as the ratio of the wavelengths 
measured at the observer and the emission (galaxy) po- 
sitions, both in their rest-frame. It can be expressed in 
terms of the tangent timelike vector to the observer (la- 
belled O) and emitting galaxy (labelled E), m m , and the 
tangent vector k^ 1 to the null geodesic joining E to O as 



1 + z 



(Al) 



We want to express the redshift and its time drift at first 
order in the perturbations around a Friedmann-Lemaitre 
spacetime with general metric ds 2 = g^dx^dx" = 
a 2 g^ v dx^dx v with 



g tiv dx^dx v 



: -(1 + 2$)dry 2 + + h ij )dx i dx j , (A2) 

2tyjij, i.e. we chose to work in the Newto- 
nian gauge and have neglected the effect of gravity waves. 

It is clear that if k^ is the tangent vector to a null 
geodesic of g^ v then fc M = a 2 k^ is the tangent vector to 



where hij — 



a null geodesic of g^ v . Decomposing fc M as k^ — E(l + 
M, e l + Se l ) k» = Ea~ 2 (l + M, e l + Se l ), where E 
is a constant, the geodesic equation reduces to (see e.g. 
Ref pj] for details) 



AM 
~d\ 



(A3) 



where A is an affine parameter along the null geodesic, 
dM/dX = M' + e l diM and a prime refers to a derivative 
with respect to the conformal time. 

At first order in the perturbations, the vector field u M 
is explicitly given by = a(— 1 — fj). We then deduce 
that fc"u M = E[-l-M-$ + e l v l ] /a. Thus 



(1 + *) 



{l-[M 



e v 



(A4) 



Integrating Eq. (|A3|) . one derives that [M] 



-2[4], 



/ _ kKj eiej ) dA. Therefore, the redshift writes 



J 



(1 + *) 



a(r) c 



1 + [$ + e^j]. 



($' + $') [a;(ry),?7]d7? 



{i+m:}. 



(A5) 



where we have shifted to the conformal time. This equa- 
tion indeed mimics exactly the standard Sachs- Wolfe for- 
mula [3. 

At the background level, the observer and emitter are 
comoving so that their proper time corresponds to the 
cosmic time. It follows that 8rj E = Sr]o = St a /ao, so 
Eq. (|A5|) implies Sz = (1 + z)(TLq — H B )5r]o/a>o, where 
H = a' I a. Shifting back to cosmic time, we get the 
standard expression for the time drift of a source located 
at redshift z, that is Eq. (TTJ). 

At first order, one has to take into account the motion 
of the observer and the emitter, as well as the metric per- 
turbations. This will manifest in the difference between 
the cosmic time and proper time. 

If at the proper time tq, the observer was located in 
(x , 770) and had a proper velocity v Q , then at a proper 
time To + St 

1. the cosmic time is St — (1 + $ )Sto = ao(l + 
<!> )<5?7o, up to terms in v 2 /c 2 and 

2. the observer has moved to x' = x Q + v Sijo so 



that he is located in 

1 - $ 

Oo'^o) = Oo^o) + (« ,1) —St 

a Q 



(A6) 



The null geodesic is still evaluated at the background 
level so that x(rj) = x' + e(r)' — rj), where e is the di- 
rection of observation, and the emitter is now located 
x' = x E + v E Sr] E . This implies that 



SVm = l + e ' VE 5 Va ~ [1 + e.{v E - v o )]St] . (A7) 
I + e.v„ 



Thus, plugging these new positions in Eq. (|A5[) . we ob- 
tain that z = z(r]o, z) + ((x ,rjo,e; z) with 



((x ,r)o,e;z) = -<S> o z(i] , z) + (1 + z) 



e.v - * 



(A3) 



where a dot refers to a derivative with respect to observer 
proper time. This expression gives the full redshift drift 
at first order in the metric perturbations and in v/c. 



